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' Abstract. The axisymmetric component of the large-scale solar magnetic fields has a pro- 

nounced poleward branch at higher latitudes. In order to clarify the origin of this branch 
we construct an axisymmetric model of the passive transport of the mean poloidal mag- 
netic field in the convective zone, including meridional circulation, anisotropic diffusivity, 
ON ' turbulent pumping and density pumping. For realistic values of the transport coefficients 

we find that diffusivity is prevalent, and the latitudinal distribution of the field at the 
surface simply reflects the conditions at the bottom of the convective zone. Pumping 
effects concentrate the field to the bottom of the convective zone; a significant part of this 
pumping occurs in a shallow subsurface layer, normally not resolved in dynamo models. 
The phase delay of the surface poloidal field relative to the bottom poloidal field is found 
■ to be small. These results support the double dynamo wave models, may be compatible 

with some form of a mixed transport scenario, and exclude the passive transport theory 
for the origin of the polar branch. 

oo 

ON 

• 1. Introduction 

The time-latitude distribution of the axisymmetric component of the large- 



scale solar magnetic fields is known to have a pronounced poleward branch 



i at higher latitudes (Stenflo, 1988, 1991, 1994; Stenflo and Giidel, 1988; Ribes 

and Bonnefond, 1990; Mouradian and Soru-Escaut, 1991). This branch is 
also present in the butterfly diagram of a number of tracers of the magnetic 
field such as quiescent filaments, polar faculae, or the coronal green line 



(Callebaut and Makarov, 1992; Makarov and Sivaraman, 1989; Leroy and 
Noens, 1983). 

Several conflicting explanations exist for this polar branch. Noting that 
the separation latitude of the two branches, 30-40°, approximately coincides 
with the latitude where the radial differential rotation changes sign accord- 
ing to helioseismology, one group of theories interprets it as the surface 
reflection of a high-latitude poleward propagating dynamo wave, coexisting 
with the low-latitude equatorward wave (Gilman, Morrow, and De Luca, 
1989; Belvedere, Pidatella, and Proctor, 1990, 1991). The Parker- Yoshimura 
rule of sign (Belvedere, 1985) then naturally leads to the correct directions 
of propagation if a is negative, as expected in the lower overshooting layer 
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where the dynamo should operate. (In fact this is only so for certain lati- 
tudinal profiles of a, cf. Schmitt, 1993.) The most modern version of these 
models, also incorporating some transport effects, is due to Riidiger and 
Brandenburg (1995). In what follows we will refer to these models as double 
wave models. 

An alternative approach to the problem of the origin of the poleward 
drift is to interpret it in terms of magnetic transport processes in the solar 
photosphere and convective zone. The common ancestor of all such mod- 
els was the now classic Babcock-Leighton interpretation of the solar cycle 
(Babcock, 1961, Leighton, 1964). The Babcock-Leighton model was what 
one may call a mixed transport model in as much as the poloidal fields were 
brought to the surface in a concentrated form in active regions, and there- 
after they were passively transported to the poles by transport processes 
(diffusion and meridional circulation). In the more recent mixed transport 
models (Wang, Sheeley, and Nash, 1991; Choudhuri, Schiissler, and Dik- 
pati, 1995) meridional circulation plays the main role in transporting the 
weak fields to the poles, and is also supposed to be responsible for the 
equatorward drift of toroidal fields in the dynamo layer. (In the original 
Babcock-Leighton theory the equatorward branch was still supposed to be 
due to a dynamo wave.) 

A third type of models may be called passive transport models. These 
are a variation of the mixed transport models where the "active source" 
of the weak fields in the form of active regions is also substituted by a 
passive transport mechanism. Dikpati & Choudhuri (1994, 1995) suggested 
such a scenario wherein the poloidal fields are brought to the surface at low 
latitudes by meridional circulation. In these models the origin of the equa- 
torward drift is not specified; instead, it is just given as a lower boundary 
condition. The resulting butterfly diagrams for the poloidal field were in 
good agreement with the observations provided that the value of the turbu- 
lent magnetic diffusivity used in the calculations was artifcially reduced to 
about 10 km 2 /s, i.e. nearly 2 orders of magnitude lower than expected. 

In the present paper we generalize the models of Dikpati & Choudhuri 
(1994, 1995) to include all turbulent transport effects such as anisotropic 
turbulent diffusion, turbulent pumping, and density pumping, as well as 
the meridional circulation. We construct models with realistic values of the 
transport coefficients in order to see whether the horizontal component of 
pumping can support the meridional circulation in producing a poleward 
migration despite the strong diffusive link between the surface and the bot- 
tom of the convective zone. (Such a possibility was suggested by Krivodub- 
skij and Kichatinov, 1991 and Kichatinov, 1993.) Another motivation for 
this study is to check to what extent is the radial pumping able to con- 
centrate the fields to the bottom of the convective zone, thereby restricting 
dynamo action to that layer, and to reduce the surface field strength to the 
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observed values, as suggested by Schiissler (1984), Petrovay and Szakaly 
(1993), and Petrovay (1994a). 



2. The model 

We calculate the mean poloidal magnetic field in the solar convective zone 
as a function of space and time, assuming axial symmetry and antisymmetry 
to the equatorial plane. No dynamo is supposed to operate in the convective 
zone, i.e. the a-effect is neglected. The dynamo layer below is not included 
in our computational volume; instead, the processes operating there are 
assumed to influence our model via the lower boundary conditions. 

The equations of passive magnetic field transport with these assumptions 
are (Petrovay, 1994b): 

d t A^ = aeed^A^ + a er 3 e 9 r A^ + a rr ^^ + 

aedgA^ + a r d r A^, + a A^ (1) 

(B e ) = -A+fr - d r A$ (B r ) = d e A^/r + cot OA^/r (2) 

where 

aee = he = Pee/r 2 a 9r = b dr = 2(3 6r /r 

ae = (70 + le~ U e )/r + {pee cot 6 - p er )/r 2 (4) 
a r = (7r + It - U r ) + (Prr + Pee + Per cot 0)/r (5) 



ao = [(7r + 7r - U r ) + cot 0(70 + 7e - U e )]/r 

-[Prr + Pee cot 2 6 + Pe r (l + cot 0)]/r 2 . (6) 

B is the mean magnetic flux density, A^ is the azimuthal component of the 
mean vector potential, U is the velocity of meridional circulation, is the 
anisotropic turbulent diffusivity, and 7 and 7 are the normal and anomalous 
components of the pumping (incorporating both the density gradient and 
the gradient of turbulence intensity, Moffatt, 1983). Note that in another 
widespread notation 7 and 7 are expressed by components of the anti- 
symmetric and symmetric parts of the a-tensor. For U we use the same 
expression as given in Dikpati and Choudhuri (1994). For the diffusivity 
and pumping standard mean-field theory expressions are used; the formulae 
are collected in the Appendix. For the evaluation of those expressions we 
use the UKX convective zone model (Unno, Kondo, and Xiong, 1985). The 
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differential rotation is assumed to be independent of depth in the computa- 
tional regime. 

For the solution of Eq. (|l|) we employ the following boundary conditions. 
At the bottom of the computational regime A<f> = A ( f ) (8,t) was explicitly 
given, as mentioned above. For the other boundaries we set 

A<f, = at 9 = (7) 

d e A^ =0 at 9 = vr/2 (8) 

d r A^ = —Afyjr at the surface (9) 

The latter condition implies that the field is vertical at the surface (Yoshimu- 
ra, 1975; Brandenburg, 1994). Note that these simple boundary conditions 
are not necessarily unrealistic: the field is intermittent, and if most of the 
flux in the photosphere is present in flux tubes exceeding 10 18 Mx then 
Eq. @ should apply owing to the strong buoyancy of these tubes. If, on the 
other hand, the flux is mainly present in thinner tubes, then a potential field 
boundary condition would be better (cf. Petrovay, 1994b for a discussion 
of this problem). We performed some test runs with different forms of the 
boundary condition to find that this choice does not have a great influence 
on the results. 

Equation (|l|) was solved by the alternating direction implicit (ADI) 
method on a 64x64 grid. The grid we use is generally uniform in both and 
r. However, immediately below the surface the scale heights of the density 
and of the turbulent velocity are very small. A grid uniform in r is unable 
to resolve this fine structure and the corresponding potentially important 
pumping effects. For this reason, in some calculations a non-uniform grid is 
employed, resolving these layers. (In practice this is realized by introducing 
a new variable f = c\ tanh(c2r — C3), rewriting eq. ([!]) in f, and solving it 
on a grid uniform in f.) 



3. Results 

3.1. Stationary dipole 

It is instructive to consider the case when the lower boundary condition cor- 
responds to the vector potential at fixed r = 500 Mm of a time-independent 
dipole field centered on the solar center. The cross section of the resulting 
field configuration after reaching a stationary state is shown in Fig. [j]. Field 
lines are nearly radial, as a consequence of diffusivity. The effects of rotation 
and meridional circulation are negligible (Fig. |2|). The overall field strength 
varies only by a factor of 2 between the surface and the bottom of the con- 
vective zone. The horizontal field component, however, increases by three 
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Figure 1. Field line configuration in a meridional plane for the stationary dipole model 



orders of magnitude from the surface to the bottom, as a consequence of 
the pumping (Fig. ||). This is the only case where the effect of rotation is 
not quite negligible; besides, it is apparent that a proper resolution of the 
shallow layers reduces the surface value by nearly one order of magnitude. 
The overall field structure is however not greatly distorted by a uniform 
grid, neither here, nor in any of the other cases studied. 



3.2. Oscillating dipole 



The case when the lower boundary condition is a dipole with a dipole 
moment oscillating sinusoidally around zero was also considered, in order 
to see whether transport effects can turn this "standing dynamo wave" into 
a migrating wave. The result was negative: owing to the overwhelming dif- 
fusivity, the surface field pattern simply reflects the pattern at the bottom 
with a slight time delay. This finding is contrary to the result of Kichatinov 
(1993) who in a simplified model found that the effect of pumping, while 
weak, is sufficient to turn a high-latitude standing wave component into a 
migrating wave. While differences in the details of the models may also con- 
tribute to the different results, we believe that perhaps the most important 
difference is that in Kichatinov's model the region of flux transport spa- 
tially coincides with the region of the dynamo (i.e. it is a convection zone 
dynamo) . 
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Figure 2. Strength of the radial field component as a function of latitude near the surface 
for the stationary dipole model (solid). Dashed: Same at lower boundary. Dots: Model 
with rotation and circulation switched off. 

3.3. Migrating field 

The lower boundary condition here was 

A* = F(6, t; Ai,ri, Ai, 0) + F(9, t; A 2 , T 2 , A 2 , S 2 ) (10) 
where 

F(e,t;A u T l ,Xi,S l ) = Aiil + expiTiiTi/A-e)}}- 1 x 

cos[ujt + 27r/X i (Tr/2-e) + 5 l ] (11) 

with oj = 9.1 ■ 10~ 9 /s the cycle frequency. 

The butterfly diagrams at the surface essentially reflect those at the 
bottom of the convective zone (Figs. || and ||). A high-latitude poleward 
branch and a low-latitude equatorward branch can thus only be produced 
at the surface if such a pattern is given as input at the bottom. Figure ^ 
shows that the field lines do not show a complicated topology inside the 
convective zone, although they are not exactly vertical either as the time 
dependence does not allow sufficient time for diffusion to smooth them. 
As a consequence, the strong downwards pumping of the horizontal field 
component is to some extent also "felt" by the vertical component, and, 
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Figure 3. Strength of the horizontal field component as a function of depth at = 45° 
for the stationary dipole model (solid). Dots: Same with a uniform grid (shallow layers 
not resolved). Dashed: Model with rotation and circulation switched off. 

in contrast to the stationary dipole case, the overall field strength is also 
significantly concentrated to the bottom of the convective zone (Fig. 0) . 

4. Conclusion 

In our axisymmetric model of the passive transport of the large-scale mean 
poloidal magnetic field in the solar convective zone we found that the lat- 
itudinal distribution of the field at the surface reflects the conditions at 
the bottom of the convective zone, i.e. in this regard the convective zone 
behaves as a "steamy window". This is due to the fact that with realis- 
tic values of the transport coefficients diffusivity is prevalent over all other 
effects. Passive transport theories of the origin of the poleward branch of 
the solar butterfly diagram are thus not viable, while double wave mod- 
els are supported by these results. Owing to the large diffusivity the phase 
delay of the surface poloidal field relative to the bottom poloidal field is 
minimal. Coupled with the approximately tt phase difference between the 
toroidal and poloidal field components for negative a value (Dikpati and 
Choudhuri, 1995), this could bring the double wave models to accordance 
with the observed phase relationship. 
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Figure 4- Butterfly diagram (i.e. B r vs. 9 and t) near the surface {top) and at the 
bottom of the convective zone {bottom) for a model with an equatorward propagating 
wave (ri = 0, Ai = tt/2, A 2 = 0). 



A mixed transport scenario of the kind proposed by Wang, Sheeley, and 
Nash (1991) is not ruled out (note that our Figs. |l| and ||| agree with the 
field structure suggested in that paper) , but the large diffusivity implies that 
the photospheric patterns should pervade the whole convective zone, being 
continuously reprocessed through it. Note that in fact the active source used 
in those models may also extend to deeper layers of the zone, as emerging 
magnetic loops are thought to shred a significant amount of flux during their 
rise (Petrovay and Szakaly, 1993; Moreno-Insertis, Caligari, and Schiissler, 
1995; Petrovay and Moreno-Insertis, 1997). 

The vertical distribution of the field (Figs. || and |?]) is strongly concen- 
trated to the bottom of the convective zone owing to pumping effects. This 
is in agreement with the earlier proposal of Schiissler (1984) and the one- 
dimensional results of Petrovay and Szakaly (1993) and Petrovay (1994a). 
The effective pumping may contribute to restricting the dynamo to the 
overshooting layer and it can reduce the surface field by about an order of 
magnitude, thereby offering a solution to the overtly high surface poloidal 
field values found by Riidiger and Brandenburg (1995) in their double wave 
model. Note that a significant part of this pumping effect occurs in a shallow 
subsurface layer, normally not resolved in dynamo models. 
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Figure 5. Butterfly diagram near the surface (top) and at the bottom of the convective 
zone (bottom) for a migrating double wave model (A1/A2 = 4, Ti = — T2 = 8, Ai = 
-A 2 = 4tt/9, 61 = 0, 82 = tt/4). 

In summary, our results support the double wave models, may be com- 
patible with some form of the mixed transport scenario, and exclude the 
passive transport theory. 

In order to decide between the double wave and mixed transport mod- 
els, two key issues should be solved in the future. First, note that in the 
mixed transport model the direction of migration is opposite for poloidal 
and toroidal fields but independent of latitude. Thus, if e.g. a clear signa- 
ture of migrating high-latitude toroidal fields were found, this could solve 
the problem in either way, depending on the direction of migration. Harvey 
(1994) claims that high latitude ephemeral active regions show an equator- 
ward migration, while Callebaut and Makarov (1992) claim that at least 
50 % of polar faculae (well known for their poleward drift) correspond to 
dipoles with a preferential east-west orientation, thus forming part of the 
toroidal field. It has even been claimed that the highest-latitude part of the 
sunspot butterfly diagram also shows a poleward drift (Becker, 1959). A 
clarification of this issue would clearly be important. 

The role of causality relations between the two branches should also be 
clarified in both theories. Surges in low-latitude activity are followed by 
poleward surges of the high-latitude field, thus supporting the mixed trans- 
port scenario. On the other hand, the well-known fact that the level of 
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Figure 6. Field line configurations in a meridional plane at 4 equidistant cycle phases for 
the migrating double wave model 



low-latitude activity in a sunspot cycle can be predicted from high-latitude 
fields at the end of the previous cycle suggests a causal relationship of the 
opposite sense. It is not impossible that such causal relations can be accom- 
modated in both scenarios, but this remains to be seen. 
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Figure 7. Total field strength as a function of depth at 9 = 45° at different cycle phases 
for a migrating double wave model 

Appendix 

Expressions for the transport coefficients 

The UKX model assumed fixed anisotropy P rr /Pee = 2 for the diffusivity 
tensor. We assume that the effect of rotation on the diffusivity remains the 
same as for quasi-isotropic turbulence (Kichatinov, 1988): 

Pee = 3/3(03 + sin 2 9) P H = 3/?0 3 
Prr = 3/?(2<fe + 2 cos 2 9) for = ~ \Pfo sin(20) 



Pe<p = P^r = (12) 

1 /ft 2 + 3 \ 

M^*) = ^[^ arctan 0,-31 (13) 

i (a 2 -i \ 

^ n *^ = Wi\n — arctana * + 1 ) ( 14 ) 

P = \x 2 t (15) 

fl, = tQ/tt (16) 

r = Max {Hp/x; 1000 s} (17) 
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where x is the r.m.s. radial component of the turbulent velocity in the UKX 
model, Hp is the pressure scale height, and f2 = £1(9) is the surface plas- 
ma differential rotation law. Note that the expression of (3 used is valid 
for non-cellular flows only (free random walk). The cellular nature of the 
supergranular flow on the solar surface may reduce the diffusivity there 
(Ruzmaikin and Molchanov, 1997), but even this reduced value is well in 
excess of 100km 2 /s and the reduction should only be confined to the pho- 
tosphere, as a similar cellular flow is not expected in the deeper layers. 

The general expressions of the pumping velocities for weak magnetic 
fields are 



le + le= ( k p^^ + k p) Pordrp/p 



H 2 

+K v ^- l (d e l3 ee + 2(3 9r )/r + K v d r 



(18) 



( 



K, 



3 Al 



K p ) f3 rr d r p/p + Kv-^g^drfi, 



+K V (def3 er + r 
+Khn ^Aj+T x 



- M jr 



(19) 



Al 



3^2 



l + {4>2l4>z)8in 2 6 



1/2 



1 + 



1 + {falfa) sin 2 



'1/2 



(20) 
(21) 



where K v = 0.6, K p = 0.15, Kp = 0.55 and K im = 0.1 account for (small) 
higher-order and intermittency corrections (Petrovay, 1994b, Petrovay and 
Zsargo, 1998); these corrections are of minor importance. Note that the 
magnetic field strength does not exceed about 10 G at any point in our 
computational volume (if we assume that the surface fields are about 1 G 
as observed); so buoyancy effects were altogether neglected here. 
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